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ABSTRACT 

Deep Chandra ACIS observations of the region around the putative pulsar, CXOU J061705.3+222127, 
in the supernova remnant IC443 reveal an ~5"-radius ring-like structure surrounding the pulsar and 
a jet-like feature oriented roughly north-south across the ring and through the pulsar’s location at 
06 h 17 m 5.200 s +22°21 , 27.52" (J2000.0 coordinates). The observations further confirm that (1) the spec¬ 
trum and flux of the central object are consistent with a rotation-powered pulsar, (2) the non-thermal 
spectrum and morphology of the surrounding nebula are consistent with a pulsar wind and, (3) the 
spectrum at greater distances is consistent with thermal emission from the supernova remnant. The 
cometary shape of the nebula, suggesting motion towards the southwest, appears to be subsonic: There 
is no evidence either spectrally or morphologically for a bow shock or contact discontinuity; the nearly 
circular ring is not distorted by motion through the ambient medium; and the shape near the apex of 
the nebula is narrow. Comparing this observation with previous observations of the same target, we 
set a 99% confidence upper limit to the proper motion of CXOU J061705.3+222127 to be less than 
44 mas yr~ 4 (310 km s _1 for a distance of 1.5 kpc), with the best-fit (but not statistically significant) 
projected direction toward the west. 

Subject headings: ISM: individual (G189.22+2.90, IC443) — X-rays: individual 

(CXOU J061705.3+222127) — stars: neutron — supernova remnants — X-rays: ISM 


1. INTRODUCTION 

Pulsars lose their rotational energy in the form of mag¬ 
netized relativistic winds. Being confined by their environ¬ 
ment, these winds produce luminous pulsar wind nebulae 
(PWNe) with shapes and spectra determined by the pul¬ 
sar emission geometry, magnetic field, and particle energy 
distribution as well as by the properties of the ambient 
medium and th e motion of the pulsar relat ive to its sur¬ 
roundings (e.s.. iKargaltsev fe Pavlovl 12008k . This is par¬ 
ticularly true of pulsars still within their natal supernova 
remnant (SNR). PWNe evolve as they first expand into 
unshocked SN ejecta, then interact with the SNR reverse 
shock, the hot SN debris, and finally the surrounding inter¬ 
stella r medium (see, e.g., the review bv lGaensler fe Slanel . 
l2006l l. Moreover, the pulsar may obtain a high space veloc¬ 
ity at birth, move supersonically through its surroundings, 
and produce a cometary PWN with a bow sh ock that may 
be vi sible in X-rays and optical emission lines dBucciantinil . 
[ 20021 ) . 

A particularly intriguing example of the PWN phe¬ 
nomena is CXOU J061705.3+222127 (also known as 
G189.22+2.90 and as 1SAX J0617.1+2221; hereafter 


,10617) located at the southern edge of the SNR shell 
IC44 3 where the she ll interacts with a molecular cloud 
('e.g. dLee et al.l . l2008lf . This alignment suggests this PWN 
is entering a rare transitional stage where the overall mor¬ 
phology of the PWN is affected by its interaction with 
the SNR shell. While it shows a dist inct cometary struc¬ 
ture suggesting supersonic motion (lOlbert et all 12001k 
iBocchino fe Bvko1l200ll ; iGaensler et al.l . 120061) . the PWN 


is oriented ~50° away from the direction expected if the 
pulsar originated at the center of the IC 443 remnant. This 
casts some doubt as to the physic al association of this o b- 
ject and IC443 as pointed out bv lGaensler et alJ (|2006h . 

A hard, non-thermal, X-ray source locat ed in the south¬ 
ern po rtion of IC 443 was first discovered bv lKeohane et all 
(|l997l) with ASCA. The source was l ater observed 
with Be ppoSAX ( Bocchino fe Bvkovl . [2000T) . with XMM- 
Newton (|Bocchino fe Bvkovl 200lt) . and with Chandra for 
~10 ks using t he Advanced CCD Imaging Spectrometer 
(ACIS) I-array (lOlbert et ah . 20071) and again for ^37.5 ks 
using the S-array ( Gaensler et all 20061) . 

The latter three observations had sufficient angular res- 
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olution to study the pulsaiH and many details of the PWN 
and surroundings. All three studies noted the soft ther¬ 
mal spectrum of the central point source is co nsistent with 
emissi on from the surface of a neutron star. lOlbert et al.l 
(l200ll) first identified the comet-shaped morphology of the 
nebula and argued it could be explained by supersonic mo¬ 
tion causing a pulsar wind to terminate in a bow shock 
and flow downstream in a synchrotron tail. The deeper 
observation obtained bv IGaensler et al.l (2006;) revealed a 
surface brightness plateau extending to ~25" downstream 
of J0617 which they interpreted as the backward termina¬ 
tion shock. They used the locations of the forward and 
backward shocks to deduce a low Mach number, ~1.2, for 
the flow and argued this was indicative of motion throug h 
relativ ely hot (shocked) supernova debris. IGaensler et, al.l 
(1200611 also identified two compact kno ts of emission a few 
arcsec onds north and south of J0617. iBocchino fe Bvkovl 
(l200ll) showed that the hard power-law spectrum at the 
core of the nebula softens with distance from J0617 and is 
consistent with synchrotron cooling models. 

Here, analysis of a much deeper, ^152 ks, ACIS S-array 
observation is presented. An overview of the observation 
and data reduction methods is given in Section [5] followed 
by analysis of the X-ray morphology, a comparison of the 
X-ray to optical emission line images, and analysis of the 
X-ray (and radio) spectrum at high spatial resolution (Sec¬ 
tion [3]). The bright source, J0617, is confirmed to display 
the physical properties of a pulsar (though no pulsations 
are detected) that powers a bright, cometary-like, roughly 
2'xl.5 / PWN. Within the nebula is a nearly circular, ~5" 
radius, ring of enhanced X-ray emission enclosing a jet-like 
feature oriented roughly north-south and passing through 
the position of J0617. The results are then interpreted 
(Section |5j) in terms of synchrotron emission from rela¬ 
tivistic particles generated by the pulsar and confined by 
the surrounding medium. 

A distance of d = 1.5 kpc (l'=0 . 42 pc ) to J0617 is 
adopted following IWelsh fe Sallmanl (120031) but note, as 
pointed out above, there remains some doubt as to the 
physical association of J0617 and the SNR IC443. There¬ 
fore, distances are given in units of di .5 = d/ 1.5 kpc where 
appropriate. 


2. DATA & METHODS 
2.1. X-ray Observations 

J0617 was observed with the Chandra X-ray Observa¬ 
tory Advanced CCD Imaging Spectrometer on 2012 Febru¬ 
ary 6 (ObsID 13736) for ~107.6 ks and again on 2012 
February 8 (ObsID 14385) for 44.5 ks. Both observations 
used an identical instrument configuration which placed 
J0617 near the back-illuminated CCD (S3) aimpoint. The 
observations were taken in full-frame timed exposure mode 
using the standard 3.2 s frame time and the VFAINT 
telemetry format. 

The data were reprocessed using the Chandra X-ray 
Center script acis_process_events (CIAO v. 4.5) to 
apply the latest time-dependent gain correction (CalDB 
4.5.5) and to flag potential cosmic-ray background events. 
Only events imaged on the back-illuminated CCD S3 are 


used in this work. Point-like X-ray sources in the field were 
identified using the sour ce-findin g util ity in the lextrct 
image-analysis program dTennant . 120061 ). The locations of 
sources identified in the two event lists were compared and 
were found to be in sufficient agreement for most purposes 
so that the two event lists were simply combined using the 
FTOOLS utility fmerge to produce a 152 ks duration fi¬ 
lial event list. In addition, a sub-pixel image was created 
by rebinning the events to 1/2 of the nominal (0.492") 
ACIS pixel size in order to examine the finer details of 
the X-ray morphology in the immediate vicinity of J0617 
(section EH) . 

Spectral fitting reported here uses t he X SPEC, 
v. 11.3.2g, spectral-fitting package (lArnau dL Il996l l. Spec¬ 
tral redistribution matrices and ancillary response func¬ 
tions were generated using the CIAO tools mkacisrmf, 
mkarf and mkwarf. Spectral analysis is confined to events 
within the 0.5—8.0 keV energy range. The energy range 
spans only 0.5—5.0 keV in those cases where a soft thermal- 
dominated spectrum produces a negligible number of high- 
energy source counts as is the case outside the PWN re¬ 
gion. 

The X-ray surface brightness and X-ray color or hard¬ 
ness ratio vary noticeably across the ~ 8 'x 8 ' S3 h eld o f 
view. The ‘contour binning’ method of ISandersI ( 2006i ) 
was therefore applied to the merged data set in order to 
map the physical properties, deduced from spectral model 
fitting, in the regions containing J0617, its PWN, and 
portions of the surrounding diffuse thermal plasma. The 
contour binning method defines spatial bins (for spectral 
analysis) which cover contiguous regions of similar sur¬ 
face brightness and is motivated by the fact that surface 
brightness variations often follow changes in physical prop¬ 
erties, such as those at shock fronts, in the X-ray-emitting 
medium. 

For this purpose, the identified point sources were first 
removed from an image constructed from the merged event 
list and back-filled with a sampling of the adjacent back¬ 
ground using the CIAO tools roi and dmfilth following 
the diffuse emission science thread{3 The resulting image 
was then accumulatively smoothed using a signal-to-noise 
ratio threshold of 30 and contour-binned using a signal-to- 
noise threshold of 100 using the accumulate_smooth and 
contbin utilities, respectively, from ISandersI (120061) . 

An earlier observation of J0617 (ObsID 05531), taken 
2005 January 12 using a similar instrument configuration, 
is used in analysis of the proper motion of J0617 (section 

EQ . 


2.2. Optical Observations 

The analysis presented here also makes use of Ha and 
[Om] images obtained 23 Nov 2012 using the Astronomical 
Research Camera£3 E2V 230-42 CCD (~0.493" pixel" 1 ) 
mounted on the 0.9 m Southeastern Association for Re¬ 
search in Astronomy telescope at Kitt Peak Observatory. 
Seeing conditions were <;2". Contemporaneous R and 
V images were used for continuum subtraction from Ha 
and [Om], respectively. The continuum-subtracted images 
were registered to the X-ray data using 2MASS objects 


n A pulsar is inferred by the nature of its X-ray properties although no pulsations have yet been observed from the compact object. 
12 http://cxc.harvard.edu/ciao/threads/diffuse_emission/ 

13 http://www.astro-cam.com 
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common to both wavelengths. 

2.3. Radio Observations 

Archival observations of IC443 at 330 MHz, 4.8 GHz, 
and 8.4 GHz, obtained with the Very Large Array 
(VLA)E!I were also analyzed and compared to the X-ray 
data. All data were calibrated and reduced using the As¬ 
tronomical Image Processing System (AIPS) data reduc¬ 
tion packag(0. 

Multiple 330 MHz observations were carried out 
from 2005 to 2007 u sing several VLA configurations 
(| Castelletti et all 1201 ill so that the final combined image 
is sensitive to a wide range of spatial structures from ^6" 
up to ^70' and therefore samples both t he PWN and the 
surrou nding extended IC 443 SNR. See ICastelletti et al.l 
(|201ll) for further details of the 330 MHz observations 
and data reduction. The final combined image has an 
angular resolution of 17" and a background noise level of 
1.7 mjy beam -1 after correcting for primary beam at¬ 
tenuations. The average SNR contribution in the region 
surrounding J0617 is 15 mjy beam -1 . 

The 4.8 GHz observation was taken in the C configura¬ 
tion on 1997 August 26 and the 8.4 GHz observation in 
D configuration on 1997 December 31. The 4.8 GHz ob¬ 
servation was taken with two 50 MHz tunings centered at 
4.8851 GHz and 4.8351 GHz and the 8.4 GHz observation 
was taken in two 50 MHz bands centered at 8.4351 GHz 
and 8.4851 GHz. Both observations used 3C 138 as the 
flux calibrator and 0632+103 as a ph ase calibrator . These 
observations were first analyzed by lOlbert et al.l (|200lh . 
These data are only sensitive to structures on scales up 
to a little over 2' and thus are not able to sample the 
emission over the entire nebula as is possible at 330 MHz. 
Here, they were re-calibrated and reduced using standard 
data reduction techniques applied within AIPS. The data 
were self-calibrated in two rounds of phase-only calibra¬ 
tion. The final, full resolution, 4.8 GHz image rms is 
0.18 mjy beam -1 with a beam size of ~5.5"x 3.8" and 
the 8.4 GHz image rms is 0.076 mjy beam -1 with a beam 
size of ~8.1"x7.1". 

3. RESULTS 

3.1. General X-ray Morphology 

The deep Chandra image (Figure [I]) reveals J0617 as 
the brightest point-like object in the field surrounded by 
a high-surface-brightness nearly-circular ring. There are 
two ‘spokes’ apparently connecting J0617 to the ring ori¬ 
ented roughly north-south which are possibly pulsar jets. 
Exterior to the ring is a comet-shaped nebula with its ma¬ 
jor axis oriented ~50° East of North and with its apparent 
forward direction to the Southwest. The co metary nebula 
has a hard X-ray spectrum as w as noted bv lOlbert et al.l 
(j200lh and iGaensler et al.l (120061 ). Predominantly soft X- 
ray emission surrounds the roughly 2 , xl.5 / cometary neb¬ 
ula and extends throughout the S3 field of view (along 
with a few faint unrelated point-like sources). The sur¬ 
face brightness of this extended emission is spatially non- 
uniform varying by factors of 4 to 5 over spatial scales of 
order 1'. As shown by spectral analysis (section [+?!) . this 


extended diffuse emission is primarily thermal in nature 
whereas the cometary nebula, ring, and related features 
are dominated by hard non-thermal spectra. 

Analysis of these features is discussed in greater detail 
in the following subsections beginning with an estimate of 
the proper motion of J0617. 

3.1.1. Constraints on the Proper Motion of 
CXOU J061705.3+ 222127 

The current observations were taken 2582+1 days since 
the 37.5-ks observation ObsID 05531 was taken with J0617 
at the aimpoint. An earlier Chandra observation, ObsID 
00760, was too short to be useful for precise measurements 
of proper motion. Four point-like X-ray sources, in ad¬ 
dition to J0617, were identified common to both ObsID 
05531 and the deeper of the two recent observations, Ob¬ 
sID 13736. The spectra of three of these four sources have 
significant fractions of their counts above 2 keV indicat¬ 
ing a high relative probability that they are distant back¬ 
ground objects and hence fixed on the plane of the sky. 
The fourth source is softer and possibl y a foregrou n d star . 
It is listed in the USNO-Bl.O Catalog (I Monet et al.l . l2003h 
with designation 1123-0129324 with a negligible proper 
motion. Offsets (independently in R.A. and Deck) be¬ 
tween the two observations were computed by minimiz¬ 
ing the error-weighted separations among the four sources. 
After applying these relative astrometric corrections, the 
99%-conhdence contour extremum sets an upper limit to 
the proper motion <310c?i.5 km s -1 with the best-fit (but 
not statistically significant) direction toward the west. 
Specifically, the estimated motion is —22.3+33.9 mas yr -1 
in R.A. and —0.1 ± 33.9 mas yr -1 in Deck 



3 0.76 3 6.9 12 19 27 37 49 62 7 


Fig. 1.— X-ray image from the merged, 152 ks duration, Chandra- 
ACIS observation restricted to the 1.0—5.0 keV energy range and 
covering a 4.75 / x3.75 / region near the aimpoint and smoothed using 
a Gaussian kernel of 3 pixel (1.48") radius. Colors are proportional 
to the number of X-ray events per pixel using a square root scaling. 
CXOU J061705.3+222127 is the bright white object surrounded by 
a ring (light blue) and the comet-shaped PWN (red and green). 
Relatively low surface brightness structure (dark blue to black) is 
visible throughout the S3 field of view. The center of the IC 443 
supernova remnant lies ^30 / to the N of J0617. 


14 The Very Large Array of the National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

15 http://www.aips.nrao.edu 
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3.1.2. J0617, the Jet, and the Ring: Subpixel Imaging 

Figure [T] displays the 2.0—4.0 keV subpixel-sampled 
region around J0617 including the ring-like and jet-like 
structures. This energy range best distinguishes the ring 
from other features and mutes the contribution from J0617 
which is brightest at lower energies (section 13.2.21) . Fig¬ 
ure [2] shows the broad-band surface brightness radial pro¬ 
file centered on the position of J0617, extending to 12.3", 
and including only azimuthal angles between PA 200° and 
340° (±70° of West) in order to avoid contributions from 
the jet-like features. This profile was obtained from the 
1/2-pixel sampled image. Also shown is a model point 
spread function for a monochromatic source at 2.5 keV 
generated using the Chandra Ray Tracer (ChaRT) and 
MARX simulation tools. The ring appears as the peak 
at 4.8" with a (Gaussian) width of 1.2". These parame¬ 
ters were estimated by fitting the surface brightness profile 
with a combined power law plus Gaussian function. An¬ 
other peak in the radial profile appears between J0617 and 
the ring. This feature is about 1" from J0617 but adding 
an additional Gaussian component shows this feature is 
not statistically significant. 



Distance from PSR (") 

Fig. 2. — X-ray brightness radial profile centered on J0617 
and restricted to azimuthal angles on the PA range 200° to 340° 
(±70° of due W). The vertical scale is the number of counts per 
0.246” xO. 246” ACIS CCD subpixel detected in the 2.0—5.0 keV 
bandpass. The thin line is a ChaRT model of a 2.5 keV monochro¬ 
matic point source. 

The jet-like feature spans from the position of the pul¬ 
sar northward and southward until it appears to intersect 
the ring. The ring brightens at these intersections. These 
intersecti ons were described as compact emission compo¬ 
nents by iGaensler et alJ (120061) . The southern jet feature 
has a surface brightness about 75% higher than the north¬ 
ern jet in the 0.5—8.0 keV energy range. 

The center of the ring is not coincident with the location 
of J0617. A spatial model representing an annulus was fit 
to the sub-pixel X-ray image in the 2—4 keV energy range 
after ignoring data within a 2.5" radius of the center of the 
ring (due to the brightness of the pulsar) and smoothing 
the data slightly (Figure [3]). The model consists of two el¬ 
liptical Gaussians plus a constant where one Gaussian has 
a positive amplitude and the other a negative amplitude 
in order to approximate an elliptical annulus. The loca¬ 
tion of the centroid of the ring model is 06 h 17 m 05.166 s 
±22°21'28.80". This location is offset by 2.7" (0.06e?i.5 ly) 


from J0617 at PA^340°, i.e., the pulsar is roughly SSE of 
the ring center. The ring is close to circular. The best¬ 
fitting parameters result in an aspect ratio of 0.86±0.05 
and an angular eccentricity, cos -1 (6/a) = 30.7 ± 5.6 deg 
where a and 6 are the major and minor Gaussian widths. 
The PA is 84±3 deg; the ellipse major axis is nearly in the 
east-west direction. 



0.246" pixels 

Fig. 3. — X-ray subpixel image of the ring surrounding J0617 
restricted to the 2.0—4.0 keV energy range and to the region within 
ten pixels (~2.5 // ) of the center of the ring (including the pulsar) 
excluded. Solid contours indicate counts pixel -1 in the data and 
range from 1.05 to 1.55 in steps of 0.125. Dotted contours indicate 
model fits to the ring morphology using the same range of levels. 
The model is the sum of two elliptical Gaussians, one with a positive 
amplitude and one with a negative amplitude (together approximat¬ 
ing an annulus), and a constant. All parameters (centroids, widths, 
amplitudes, and position angles of the Gaussians and the constant) 
were allowed to vary in the fitting although the centroids of the two 
Gaussians were tied to the same values. Circles are shown centered 
on the best-fit centroid (denoted by the asterisk) and show that the 
overall shape of the ring is close to circular. The position of J0617 
is denoted by the red cross. 


3.1.3. Morphology of the Pulsar Wind Nebula 

Figure [4] displays the 0.5—8.0 keV brightness profile 
along the major axis (PA 50°) of the comet-shaped PWN. 
The profile has been averaged over a ±30" region perpen¬ 
dicular to this axis. The peak of the profile coincides with 
the pulsar and offsets upstream of the pulsar (toward the 
southwest) are taken as positive in this figure. The ring 
is clearly visible as secondary peaks on either side of the 
pulsar. The profile >5" upstream of the pulsar decreases 
smoothly out to at least 2' from J0617 without any dis- 
cernable structural features. As with the upstream profile, 
the downstream profile first declines steeply from the exte¬ 
rior of the ring but then it reaches a broad peak extending 
from 15" to 25" away from the pulsar. The profile then 
declines more slowly out to approximately 2! downstream 
without any additional features. 
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Fig. 4.— X-ray brightness profile along the major axis of the 
PWN. Upstream is to the right. The vertical scale is the num¬ 
ber of counts per 0.492 // x 0.492" ACIS CCD pixel detected in the 
0.5—8.0 keV bandpass in a 60" wide strip oriented perpendicular to 
and centered on the major axis. 

The X-ray image does show that the surface bright¬ 
ness distribution upstream of the pulsar is bow-shaped. 
Parabolas were fit, by visual inspection, to the surface 
brightness contours in the upstream region of a smoothed 
image (with the ordinate axis of the parabolas constrained 
to be parallel to PA 50°) as shown in Figure [5] The 
parabolas are of the form ( y/d s ) = 0.54(ir/d s ) 2 + 1 where 
d s is the distance from the location of J0617 (at the focus) 
to the vertex of the parabola. 



-4060 4080 4100 4120 4140 


Fig. 5.— X-ray image of the upstream region around the PWN. 
The image has been smoothed by three applications of a nearest- 
neighbor smooth. Contours are drawn at levels of 1.5, 3.0, 5.0, and 
10 c pixel -1 . Corresponding (color-matched) curves are parabolas 
parameterized by visual inspection to represent the contours. The 
parabolas are constrained to have their ordinates parallel to the 
apparent line of motion of J0617 at PA 50°. 

There is no evidence for a bow shock in either the Ho: 
or [Om] images of this region and there is little or no 
correlation between the X-ray emission in the upstream 
direction with the optical line emission. However, as a 
close-up image of the region shows (Figure [6]), there is Ha 
emission coincident with the jet-like X-ray feature. This is 
likely an accidental coincidence as there are many similar 


Ha features throughout the IC 443 SNR. Similarly, there 
is also a well-defined Ha wisp downstream of J0617 that 
runs roughly along the major axis of the PWN and thins 
in the downstream direction (Figure Q. This feature is 
also visible in the [O III] image. There is no X-ray fea¬ 
ture correlated with this wisp. Although its location is 
suggestive, it is unclear if this Ha wisp is related to the 
PWN as there are numerous wisps with similar morphol¬ 
ogy throughout the region that are obviously not related 
to the PWN proper. 



Fig. 6. — Detail of the ring and jet-like features in X-ray (left 
two panels) and Ha (right two panels) with X-ray contours (top 
panels) and Ha contours (bottom panels) overlaid. Panels sizes are 
l'xO.75'. 



Fig. 7.— X-ray (blue), Ha (red), and [Om] (green) color compos¬ 
ite of a 4.2'x3.4' portion of the IC443 field. The location of J0617 
is marked with a white cross. Note the Ha wisp extending east-west 
roughly along the line of apparent motion of J0617 superimposed 
on the PWN. 

3.2. X-ray Spectroscopy of Major Structures 
3.2.1. Regions Exterior to the PWN: Background SNR 

Customarily, a background region is selected from those 
regions remote from the primary source, in this case the 
PWN, in order to isolate the source spectrum for analysis. 
However, J0617 is located along the line of sight through 
the extended X-ray emission associated with the SNR 
IC 443 and this emission is highly non-uniform as noted 
above. Several regions were visually selected and analyzed 
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to quantify the spectral properties of this “background” 
throughout the field of view of the S3 chip. These regions 
are numbered 1 through 9 in Figure |8] and their spectra 
are shown in Figure [9] Assuming this emission is from 
hot plasma within the SNR, the observed variations may 
be due to a combination of different plasma temperatures, 
emission measures, elemental abundances, and interven¬ 
ing absorption. The simplest model that can account for 
all these variables in XSPEC is an absorbed optically-thin 
thermal plasma model. The spectral analysis reported in 
this section is limited to the range 0.5—5.0 keV because of 
the paucity of source photons above 5 keV. 



3800 4000 4200 4400 4600 
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Fig. 8. — Regions selected for spectral analysis of extended emis¬ 
sion. Numbering corresponds to the region numbers in Table [l] Ad¬ 
ditional small rectangular regions upstream of J0617 are analyzed in 
section ^. 2. 41 The thick irregular region around J0617 delineates the 
outermost of 7 high-surface-brightness PWN regions investigated in 
section 13.2.51 and the thin inner irregular region defines the radio- 
bright region analyzed in section 13.2.61 


It was found that the spectra of these ‘outlier’ regions 
cannot be fit satisfactorily with a single-component ab¬ 
sorbed (XSPEC model phabs) variable-abundance ther¬ 
mal emission (vmekal) model. A significant improvement 
could be made by adding either a second thermal compo¬ 
nent or a power law (powlaw) component. 

Resulting best-fit model parameters are listed in Table [I] 
for these outlier fields. These include the absorption col¬ 
umn density, nn, the power law index, T, the thermal 
component temperature, kT, the model component nor¬ 
malizations, Aft the 0.5—5.0 keV flux, /, in erg cm -2 s -1 
(uncorrected for absorption) and the x 2 fit statistic and 
number of degrees of freedom, dof. Here and elsewhere 
errors are 90% confidence extremes for a single interesting 
parameter unless otherwise noted. 



Fig. 9.— Observed spectra of the 9 numbered regions shown 
Figure [8] The spectra have been binned up to 10 channels in or¬ 
der to achieve a minimum detection significance of 10 for display 
purposes. Colors correspond to regions l:cyan, 2:gray, 3:magenta, 
4:blue, 5:black, 6:green, 7:red, 8:orange 9:dark gray as defined in 
Figure[8]and Tablc[l] Note the high flux above ~2.0 keV in the spec¬ 
trum of regions 7—9 indicating a strong hard spectral component. 
The different relative strengths of emission lines among the various 
regions are due, in part, to real differences in elemental abundances 
as the best-fitting models all have similar plasma temperatures and 
absorption columns (Table [T}. 

The first thermal component temperature is typically 
kT ~ 0.5 to 0.8 keV. It was found that the second com¬ 
ponent usually had a steep power law slope (photon index 
r ;> 2.5) or only a slightly higher temperature (kT ~ 0.8 
to 1.0 keV) suggesting that this component may represent 
the range of thermal emission temperatures expected along 
a given line of sight through the SNR rather than a true 
non-thermal component typically indicated by a power law 
shape. This is true of most of these outlier regions except 
regions 7—9 where F < 2.3 (or kT 2 > 3.5 keV). In these 
cases, it is likely that there is truly a hard non-thermal 
component even if the two -kT model is formally a better 
fit as such high temperatures are atypical of SNR spec¬ 
tra. The PWN likely contributes this non-thermal emis¬ 
sion component in regions 8 and 9 because they are imme¬ 
diately adjacent to the high-surface-brightness PWN. 

The elemental abundances were not well constrained in 
any of these models but strong emission lines of H-like O 
and Ne, and He-like blends of Ne, Mg, Si, S are present 
throughout as shown in Figure |9] whereas Fe emission is 
weak. This is indicative of a Type II core collapse SN 
origin. 

Overall, since the spatial scale of the observed brightness 
variations are of order 1', and the spectra of the different 
regions sampled also vary, none of these regions (nor any 
other region on S3) can be safely defined as a background 
to be used in fitting spectra of the PWN. Thus, the “blank 
sky” background 17 ! appropriate to this observation is the 
only background used in subsequent spectral analysis with 
the exception of the pulsar, J0617, and the jet-like and 
ring features. For these objects, the space between the 
pulsar and ring is relatively faint and local backgrounds 
were chosen from this region. 

— 14 (47rd 2 ) -1 times the emission measure, in 


16 A'pow is in units of flux density, ph keV 1 cm 2 s 1 at 1 keV and A'thm is in units of 10 
cm -3 . 

17 obtained from the CalDB and accessed following the CIAO science thread http://asc.harvard.edu/ciao/threads/acisbackground/ 
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3.2.2. Spectra and Timing of CXOU J061705.3+ 222127 

The spectrum of J0617 was extracted from a 2-pixel 
(0.984") radius circle centered at the location of the pul¬ 
sar (see also Figure ED- A nearby background was ob¬ 
tained by combining data from five 2-pixel radius circles 
located interior to the ring and excluding the jet-like fea¬ 
tures. This results in an estimated 925±35 net source 
counts in the 0.5—8.0 keV energy range. Absorbed black- 
bod y a nd ne utron star atmosphere ( Pavlov et al. . 199 5: 
IZavlinl . [2009) models (plus a power law component) were 
fit to the spectrum over the 0.5—8.0 keV energy range 
after grouping to ensure a minimum of 20 background- 
subtracted events per spectral energy bin. The spec¬ 
trum, models, and fit residuals are shown in Figure [10] 
The neutron star atmosphere model assumed a distance 
to the star of 1.5 kpc, a radius 7 ?ns = 10 km, a mass 
Mns = 1.4 Mq (corresponding to a gravitational redshift 
parameter g r = [1 — 2GMns/(-Rnsc 2 )] 1 / 2 = 0.766) and a 
magnetic field strength of 10 13 G. 




channel energy (keV) 


Fig. 10.— Observed spectrum of CXOU J061705.3+222127 ex¬ 
tracted from within an ~1 // (2-pixel) radius circle. A background, 
extracted from small regions surrounding the source, interior to the 
ring, and excluding bright regions associated with the jet-like fea¬ 
ture, has been subtracted. The spectrum has been grouped to obtain 
a minimum of 20 counts per spectral bin. The solid line denotes the 
full model (top; absorbed blackbody plus power law; bottom: ab¬ 
sorbed neutron star atmosphere plus power law), The dotted lines 
are the contributions of the power law component and the dashed 
lines are the blackbody (top) and NS A model ( bottom ) components. 
The lower panels show the respective fit residuals. 

The best-fitting model parameters, derived properties, 
and fit statistics are listed in Tabled where T°° and R°° 
are the effective temperature and radius as measured by 
a distant observer, and Lr is the 0.5—8.0 keV luminos¬ 


ity of the non-thermal (power law) component. Errors in 
Table [2] are at the lcr (68%) confidence level. Both mod¬ 
els provided acceptable fits to the data. These spectral 
properties are typical for a cooling neutron star although 
the blackbody model radius, R°° = 1.63 km, implies emis¬ 
sion from a small hotspot rather than the entire neutron 
star surface. They indicate a bolometric luminosity of 
(1—3)xl0 32 erg s” 1 . 

A search for a possible periodicity in the X-ray flux from 
J0617 was restricted to periods P > 6.5 s due to the 3.24- 
s time resolution of the ACIS instrumental setup. The 
Rayleigh test yielded a maximum power at a period of 
31. 18 s. Apply ing the appropriate probability distribution 
dGrotbl. 1197511 sets an upper limit to the power at 99% 
confidence that translates to an upper limit on the pulsed 
fraction of 36% assuming a sine-wave signal and account¬ 
ing for the background. 

3.2.3. Spectra of the Jet and Ring 

The spectrum of the southern jet-like feature was ex¬ 
tracted from a 10" x 5" region oriented N—S and centered 
at 06 h 17 m 5.25 s +22°21'25.5". The background-subtracted 
source spectrum (using the same background regions as 
for J0617) was grouped to a minimum of 10 events per 
spectral energy bin and fit using an absorbed power law. 
There were too few counts (233±21) to justify attempt¬ 
ing more complex models. The best-fitting model pa¬ 
rameters are «h = 6 .3j/ 3 '2 x 10 21 cm -2 and T = 

(x 2 /dof = 20.9/18) resulting in an absorbed 0.5—8.0 keV 
model flux of 2.l/l^'g x 10” 14 erg cm -2 s” 1 and absorption- 
corrected luminosity of 7.0)l4;o x 10 29 d 2 5 erg s” 1 . There 
were too few counts for spectral analysis of the northern 
portion of the jet. 

The spectrum of the ring was extracted from a cir¬ 
cular annulus lying between radii of 3.7" and 5.9" from 
the ring center. An absorbed power law model with pa¬ 
rameters riR = 5.2jjg;® x 10 21 cm -2 , F=1.4±0.1 produced 
an acceptable fit to the spectrum (y 2 /dof = 194.3/182). 
The resulting 0.5—8.0 keV model flux is (3.1±0.3)xl0” 13 
erg cm” 2 s” 1 and the absorption-corrected luminosity is 
l.lxl0 32 df 5 erg s” 1 . 

3.2.4. Spectra Upstream of the Pulsar Wind Nebula 

Spectra were extracted from a series of regions located 
upstream of the pulsar to search for a spectral transition 
between the non-thermal shocked PWN emission and ther¬ 
mal emission from shocked ambient medium further up¬ 
stream. Figure [8] shows the locations of these regions. 
They are 40x100 pixel 2 (0.269 arcmin 2 ) rectangles ori¬ 
ented with the long side perpendicular to the line of ap¬ 
parent motion of the pulsar. Table [3] presents the spectral 
fitting results for these regions (ordered left to right by in¬ 
creasing distance from J0617). As with the outlier regions 
discussed in section 13.2.11 absorbed two-component mod¬ 
els were applied. Both models provide similarly acceptable 
fits although the thermal plus power law is formally the 
superior model. With the exception of the region nearest 
to J0617, the temperature of the soft thermal component 
in all cases, kT ~ 0.6 — 0.7 keV, is consistent with the soft 
component deduced for the outlier regions of section O'. 2. II 
Most regions have an obvious hard component (T ^ 2.3 or 
kT 2 > 3.5 keV) that is poorly modeled by a second ther- 




































mal component. For the non-thermal power law (plus soft 
thermal) models, the trend is for the power law index to 
steepen and both the power law flux and its fraction of 
the total flux to decrease with distance from J0617. This 
is certainly true of the innermost three regions although 
the trend is less evident in the outer three regions where 
the power law component is much weaker though still re¬ 
quired in the modeling. While there is no abrupt spec¬ 
tral transition from non-thermal- to thermal-dominated 
emission, a non-thermal component is clearly present at 
least out to 2' upstream from J0617, which is consistent 
with the smoothly-varying surface brightness distribution 
of this upstream region (cf. Figure [4] and section 13.1.31) . 


flatten slightly to compensate. However, this is only a 
modest change in slope and within the errors obtained if 
the thermal model temperature were left as a free param¬ 
eter. For example, in the fourth zone, T decreases from 
1.9±0.1 to 1.8±0.1. 

Maps of several of the resulting best-fitting model 
parameters are displayed in Figure [Tl] for the single¬ 
component absorbed powerlaw model. These maps show 
that the power law slope smoothly increases with distance 
from J0617, that the intensity decreases, and that the best¬ 
fitting column density is roughly constant but decreases in 
the outermost zone. This latter trend is perhaps influenced 
by the simplicity of this single-component model. 


3.2.5. Spectra of the Pulsar Wind Nebula: Contour 

Binning 

The contour binning method of lSandcrsl (20061 1 was ap¬ 
plied to a region containing the pulsar, PWN, and im¬ 
mediate surroundings as described in jj2] By definition, 
this binning method traces the decrease in surface bright¬ 
ness with distance from J0617. Using this prescription, 
the PWN is subdivided into roughly concentric “annuli” 
with the innermost region containing J0617, the ring and 
the jet. The outermost of these regions is indicated in Fig- 
ure[8]by the thick irregular contour surrounding J0617. By 
design, each region contains roughly 10 4 X-ray events. Be¬ 
cause the extended thermal emission from the IC 443 SNR 
is so spatially non-uniform, the contour binning method 
did not satisfactorily bin the regions beyond the PWN. 

The results of model fits to the spectra of these regions 
are given in Table [4] (a 3-pixel (1.48") radius circle centered 
on J0617 was excluded from the innermost region, denoted 
Region 1 for spectral analysis; regions 2 through 7 lie se¬ 
quentially ourward). In addition to model parameters and 
derived properties listed previously (Table |T] and [3]), Ta¬ 
ble U] includes fluxes, / cor , corrected for absorption along 
the line of sight. The spectra were fit on the extended 
0.5—8.0 keV energy range and fluxes listed in Tableware 
over this same interval. 

An absorbed power law model provided a statistically- 
superior fit to the spectra of the three innermost regions. 
An absorbed two-component (powerlaw and vmekal or 
two vmekal) model was needed for the outer zones but 
one component is always hard (r < 2 or kT > 5 keV) 
and dominates the 0.5—8.0 keV flux; again, indicating a 
true non-thermal component is present as argued in sec¬ 
tion 13.2.11 For this reason, the two vmekal component 
model is not considered further. The soft thermal com¬ 
ponent in the two-component non-thermal model is gener¬ 
ally cooler, kT 0.4 keV, throughout these PWN regions 
than in the ‘outlier’ fields (Section 13.2.11) but note that 
the soft component contributes at most 5% of the flux in 
all but the outermost PWN region. The outermost region 
is similar to the outlier regions 7—9 (Table Q]) in that its 
soft thermal component temperature is 0.67 keV and its 
power-law index is ~2. Unlike these outlier regions, nearly 
all the flux in even this outermost PWN region is provided 
by the non-thermal model component. 

To test the significance of the low temperature in the 
inner zones, the fits were repeated with the vmekal model 
temperature parameter fixed to its value in the outermost 
zone, kT = 0.67 keV. This results in some additional flux 
in the 1-2 keV band which forces the power law index to 


3.2.6. Spectra of the Pulsar Wind Nebula: Comparisons 
to Radio Data 

The X-ray emission in the J0617 pulsar wind nebula is 
clearly dominated by a non-thermal component suggesting 
a synchrotron emiss i on source. Radio i mages of the region 
(jOlbert et al.l . 120011 : ICastelletti et all l201lh also show a 
cometary morphology which may imply a similar origin. 
However, the radio and X-ray data have much different 
spatial resolution and sensitivity so that a direct compar¬ 
ison is non-trivial. Specifically, the signal to noise of the 
higher frequency radio data is too low to allow for mea¬ 
surements in small regions for comparison with the X-ray 
fluxes. Instead, by cutting the upper and lower uv range 
of the 4.8 GHz and 8.4 GHz measurements to a common 
~1.7'x0.77' region (see Figure [5]) with equal coverage at 
both frequencies (and at 330 MHz), radio spectral index 
measurements of roughly the whole J0617 PWN can be ob¬ 
tained. Comparison of the coverage after cutting showed 
no apparent gaps in the coverage and thus both frequen¬ 
cies should be sensitive to emission over a similar range 
of spatial scales. The data were therefore re-imaged and 
convolved with a circular beam of size 8.5" and the images 
were blanked for emission falling below a three sigma cut¬ 
off before measurements were taken of the total flux of this 
region. This yields nearly equal fluxes at all three frequen¬ 
cies, 84.3±7.6, 87.6±3.0, and 87.0±0.6 mJy at 330 MHz, 
4.8 GHz and 8.4 GHz, respectively, and therefore a flat 
spec tral index of or = 0.01 in the radio band consi stent 
with lOlbert et al.l (1200 lh and ICastelletti et al.l (1201 it) . 

The combined 330 MHz image does not suffer from these 
sensitivity issues but does not have the high spatial reso¬ 
lution of the X-ray data. Therefore, instead of the contour 
binning regions of section l5'.2.5l a set of successively larger 
elliptical annuli approximating the innermost five contour 
binning regions was analyzed at both 330 MHz and X- 
rays. As before, the X-ray spectrum is predominantly non- 
thermal and a single absorbed powerlaw model provided 
an acceptable fit. Here, J0617 was not excluded since it 
cannot be cleanly masked out of the low-resolution radio 
image. The 330 MHz radio fluxes were calculated by in¬ 
tegrating over each region after subtracting an estimated 
contribution due to emission from the surrounding SNR 
as noted in section 12.31 Tabic [5] summarizes the geometry 
and the radio and X-ray properties for these regions. 

If the SNR-background-subtracted 330 MHz radio flux 
and the non-thermal (power-law component) X-ray emis¬ 
sion arise from the same physical mechanism, then their 
spectral shapes should be simply related. The slopes of 
these spectra are obviously different: The radio spectrum, 
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Fig. 11.— Spectral properties of the seven regions encompassing the PWN as defined by the contour binning method of ISandersI 1120061 1 
and listed in Table QT] Properties are shown for the single absorbed power law model. Region numbers listed in Tabic [4] column headings 
(1—7) are ordered from the inside outward. Each region encompasses an area of approximately uniform surface brightness. The innermost 
(region 1) is approximately 15 pixels (~7.5”) in radius and contains the ring, and the jet features (J0617 is excluded); region 2 encompasses the 
“plateau” region (see Figure [p. The panels display (clockwise from top left:) the fitted absorption column density, rrjj/10 22 cm 2 , the power 
law model photon index, F, the observed 0.5—8.0 keV intensity, ex/10 -16 erg cm -2 s —1 arcsec -2 , and the power law model normalization, 
Upow/10 4 photons kev -1 cm -2 s —1 at 1 keV. 


aR, is flat, at least on average over the entire region, while 
ax = r — 1 ~ 0.5 — 1 depending on the (elliptical) re¬ 
gion. For a power-law X-ray spectrum with a photon 
index T and normalization Kpow = 10 - 4 AT _4 photons 
(cm 2 s keV ) -1 at E = 1 keV, the energy flux spectrum is 
F Vi x oc K- 4 iz -r+1 . Extrapolation of this X-ray spectrum 
towards lower frequencies intersects the extrapolation of 
the radio spectrum, F v> r = i 7 iGHz(p/lGHz) - “ R , at fre¬ 
quency v when 

(ax - or) \ogu = log(ifr_4/EiGHz) + 8.3838a x - 4.1798, 

where v is in units of GHz and Fighz is in mJy. Insert¬ 
ing the appropriate values from Table [5j and assuming 
or = 0 everywhere (and no uncertainty in the radio flux 
densities), the estimated frequencies v are as listed in the 
last row of Table [51 In the X-ray hard (flat) inner ellip¬ 
tical regions 1 and 2 , the extrapolation crosses the radio 
spectrum at frequencies lower than 4.8 or 8.4 GHz which 
means that there should be at least two breaks in a power 
law spectrum or the radio spectrum is not flat in those 
regions. 

4. DISCUSSION 

4.1. Properties of the Central Source 

Deep Chandra /ACIS observations show the point-like 
object, CXOU J061705.3+222127, has an X-ray spec¬ 
trum dominated by a thermal component with param¬ 
eters typical of a cooling neutron star. The fit with 
a hydrogen atmosphere (NSA) model shows an effective 
temperature T^g cs 6.8 x 10 5 K and a bolometric lumi¬ 
nosity 2.6 x 10 32 erg s -1 (see Table 0. Such 

temperature and luminosity are very close to those ob¬ 
tained from t he NSA fit of the th ermal component of the 
Vela pulsar (jPavlov et all l200ll) . whose age r is about 
20-30 kyr, larger than the characteristic spin-down age 
T sc j = P/2P = 11 kyr dLvne et al.l . 119961 ). Similar tem¬ 
perature and luminosity for r ~ 30 kyr we re obtained 
in so me neutron star cooling models (e.g., iPage et all 
l2009f ). Therefore, we can expect that age of the neutron 
star in J0617 is a few times 10 4 yr. It is consistent with 
the IC443 age ~ 30 kyr estimated bv IChevalied (119991) . al¬ 


though so me authors inferred much younger ages for this 
SNR (e.g. JPetre et al.l (119881) suggested r ~ 3 kyr). 

Compared to the NSA model, the blackbody fit of the 
thermal component gives a higher temperature, Tbb ~ 
1.5 x 10 6 K, and a smaller radius, Rbb ~ 1 . 6 di .5 km, 
similar to the Vela pulsar. Such emission could emerge 
from a hot spot on the neutron star surface heated by 
precipitation of relativistic particles from the pulsar mag¬ 
netosphere, or it could reflect strong nonuniformity of the 
surface temperature due to anisotropy of heat conductiv¬ 
ity of the neutron st ar crust caused by very s trong mag¬ 
netic fields (see, e.g., iGeppcrt fe Viganbl . l2014l) . The non¬ 
detection of pulsations with P > 6.5 s does not contradict 
the small size of the emitting region because such young 
pulsars have much shorter periodJ 43 !. 

Since the NSA and blackbody models yield fits of about 
the same quality, we cannot prefer one model to the other 
based solely on the observations. However, the NSA fit 
looks somewhat more preferable because it gives reason¬ 
able cooling neutron star properties at the expected SNR 
age and distance. 

The detection of the nonthermal power-law component 
in the neutron star spectrum and the presence of the non¬ 
thermal X-ray nebula around the neutron star show that 
CXOU J061705.3+222127 is a rotation-powered pulsar. 
The nonthermal component, presumably emitted from the 
pulsar’s magnetosphere, dominates the pulsar spectrum at 
E ^ 1.7 keV. The spectral slope and the “isotropic lumi¬ 
nosity” of the nonthermal component (Lp = 47 t d 2 /pi° r ) are 
also close to those of the Vela pulsaiFl. for both thermal 
component models used in fit. 

Since pulsations have not been detected, we do not 
know the pulsar’s spin-down properties, E and r s( j. We, 
however, can constrain them from a comparison of the 
PWN 0.5-8 keV luminosity L pwn ~ 1.4 x 10 33 d 2 5 erg 
s -1 (given by the sum of the non-thermal components 
of the 7 contour-binning regions of Section 13. 2. 51) with 
the X-ray luminosities of PWNe produced by pulsars with 
kno wn E and t s h- Assuming di n = 1 and using Figure 5 
from iKarealtsev fe Pavlov! (120081). Figure 2 from lLi et al.l 
(|2008h . and Figure 3 from lKargaltsev et al.l (|2012f ). we ob¬ 
tain a range of possible spin-down luminosities E k, (1- 


115 Periods larger than 6.5 s are observed in magnetars, but this neutron star does not show any magnetar properties. 

12 Note a typo in IPavlov et al.l (120011 1: The power-law luminosities in Table 1 should be in units of 10 31 (not 10 32 ) erg s 1 . 
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30) x 10 36 erg s , with the most likely value around 
3 x 10 36 erg s -1 . Similarly, using Figure 5 from iLi et al.l 
iH), we obtain a range of possible spin-down ages, 
r s d ~ 10-100 kyr, with the most likely value around 20 
kyr, consistent with the age estimates from the proper¬ 
ties of the thermal component. From these estimates we 
obtain constraints on the period, P ~ 0.1-0.6 s, and the 
surface magnetic field, B ~ (0.4-4) x 10 13 G, with most 
likely values around 0.25 s and 2 x 10 13 G, respectively. 

Note that the true pulsar’s X-ray luminosity is appar¬ 
ently much higher than the above-estimated Lr because 
the (distance-independent) ratio T pW n/Tr ~ 100-200 is 
much larger than for any k nown pulsar-PWN pair (see 
Figure 5 in iKargaltsev fe Pavlovl l2008h . This suggests 
that the pulsar beam misses the Earth, which could also 
be the explanation for non-detection of the pulsar in the 
radio and y-rays. 

J0617 is separated from the apparent center of the 
IC443 SNR by 15' in a direction nearly due south. This 
corresponds to a projected length sj_ = 2 x 10 19 di .5 cm 
and a projected velocity u_l = 640 di.^rf 1 km s -1 . 
This is below the observed upper limit (Section 13.1.11) 
of 310 di .5 km s -1 if 74 = r/10 4 yr ;> 2, which is con¬ 
sistent with the possible IC443 age and the spin-down 
age estimated above from the empirical L pwn -T sc i corre¬ 
lation. The true speed is u psr = ux(sini) -1 , where i is 
the angle between the direction of motion and the line of 
sight. Observed pulsar sp e eds r ange up to ~U000 km s -1 
(e.g., [Arzoumanian et ~aH . l2002h which constrains the an¬ 
gle: sin* ;> 0 . 647-4 dxs, assuming J0617 was born in the 
center of IC443 (e.g., i ^ 19° for 74 ~ 2, di.5 = 1). The 
angle i is larger for more typical pulsar velocities - for in¬ 
stance, 32° ^ i <; 53° for v = 400 km s -1 and a reasonable 
age range 2 <; 74 <; 3, at di .5 = 1. 

4.2. The surrounding medium 

All regions imaged by Chandra surrounding J0617 and 
its PWN have a strong thermal emission component, likely 
emitted from the IC 443 SNR. This likely means that J0617 
lies within this warm plasma, but a simple line-of-sight 
projection cannot be ruled out (and, thus, the distance 
to J0617 remains unknown). The typical temperature of 
this warm plasma is kT ~ 0.7 keV (Section 13.2.11) cor¬ 
responding to a sound speed, c s = ('ykT/ /rmn ) 1 ^ 2 ~ 
430 (fcT/0.7keV) 1 / 2 (0.6//z ) 1 / 2 km s -1 for a 7 = 5/3 gas 
and molecular weight p, ~ 0 . 6 , appropriate for fully-ionized 
medium with roughly cosmic abundances. Thus, if J0617 
is within this plasma and not just lying along the same 
line of sight, then the pulsar motion cannot be highly 
supersonic; the Mach number M = vo/c s <, 2.3 for 
Wo < 10 3 km s -1 . Here, vq is the speed with respect to 
the ambient medium; vo ~ w psr in the absence of any bulk 
motions in the plasma. Whether the speed is supersonic, 
transonic, or subsonic has a strong bearing on the shape of 
the nebula (which is also influenced by its orientation rel¬ 
ative to the line of sight) and the properties of the particle 
and electromagnetic field flow within the nebula. This is 
addressed in more detail in Section IT3l 

Assuming, for the present, that J0617 and its PWN 
are embedded in the IC443 SNR, the ambient thermal 
pressure, pthm = p(kT/pmn), where p is the ambient 
mass density, can be estimated with the aid of the vmekal 


model fit parameters. The vmekal model norm is TFthm = 
10” 14 (47r<i 2 ) _1 J n e nndV ~ 10 -14 (47ni 2 ) _1 n| f R, where 
nn is the hydrogen number density, and V is the volume 
of the emitting region. The volume can be expressed in 
terms of the area A of the region image (in arcmin 2 ) and 
the average depth l w (in units of 10 19 cm) of the region 
along the line of sight: V = 1.8 x 10 55 A^igd 2 5 cm 3 . Us¬ 
ing this expression for V and assuming that fully ionized 
hydrogen gives the main contribution, the hydrogen num¬ 
ber density and the ambient thermal pressure can thus be 
estimated as follows: 

nn = 0.39 k]/ 2 (AIi 9 )~ 1//2 cm" 3 , 
p thm ~ 2n H kT = 0.92 x lO^lT^A/ig) -172 kT dyncm -2 , 

where TFthm = 10 _4 /\_4 and kT is in units of keV. 

The nu and pthm estimates are proportional to l lg ' , 
which depend on the SNR geometry and the J0617 posi¬ 
tion within the SNR. J0617 lies, in projection, about 15' 

(2 x 10 19 di.5 c m ) from the center of the ~ 20' radius IC 443 
remnant (assuming a spherical remnant). Thus, a chord 
through the remnant along the line of sight to J0617 should 
have a length of about 2 x 10 19 cm or less, so typical ex¬ 
pected values are 1 19 ~ 1 (corresponding to ~ 7.5' in the 
sky plane). 

Using the values from spectral fits to the nine outlier 
regions, Table [T] the mean hydrogen densities (and stan¬ 
dard deviations) for the low- and high-temperature vmekal 
components are nn/lf 9 ^ = 0.5 ± 0.1 and 0.6 ± 0.2 cm -3 . 

_1 /o 

The corresponding pressures are p_g// 19 ' = 0.7 ± 0.2 for 
the low-temperature component (where pthm = 10 -9 p_g 
dyn cm -2 ). For the high-temperature component in re¬ 
gions 1 through 6 the pressures are 1.4 ± 0.4, but they 
average 5 ± 2 for regions 7 through 9 because of the 
much higher temperatures determined for those regions, 
averaging 4 keV compared to only 0.8 keV for the high- 
temperature component in regions 1—6 (and 0.6 keV for 
the low-temperature component in all regions). 

4.3. Morphology of the Ring and Cometary Nebula 

The cometary shape observed from many PWNe is 
generally attributed to motion of the system relative to 
the ambient medium. In the case of powerful, young, 
subsonically-moving pulsars such as the Crab, a distinct 
jet-torus structure may be present near the pulsar related 
to axial symmetry in the wind caused by rotation of the 
pulsar. In other systems, particularly older pulsars mov¬ 
ing supersonically through a cold interstellar medium, the 
PWN structure consists of a forward bowshock and an ex¬ 
tended tail. 

J0617 appears to have some attributes of both these ex¬ 
tremes. It has a distinct, nearly circular, ^"-radius cir- 
cumstellar ring encircling jet-like features and a ~2'xl.5 / 
comet-shaped nebula. In analogy to the Crab, the 
ring may define a termination shock (TS) where the ul- 
trarelativistic free-flowing equatorial wind slows to sub- 
rclati vistic speeds, referred to as the inner ring in the 
Crab (jWeisskopf et all I2000T) . or the ring may be a torus 
of shocked pulsar wind similar to the one that lies some 
distance beyond the inner ring of the Crab. Since in both 
cases the symmetry axis of the J0617 ring should be the 
spin axis of the pulsar, the ellipticity of the ring implies 
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the angle £ between the line of sight and the spin axis 
is about 30° (see Section 13.1.21) . If we assume that the 
spin axis coincides wit h the direction of mo tion, similar 
to many other pulsars (I Johnston et al.l . [20051) . and J0617 
originated at the center of IC 443, then the pulsar velocity 
can be estimated as u psr ~ 400-600 km s _1 , for di. 5 = 1 
and T4 ~ 2-3 (see Section |4~T1) . comparable to the speed 
of sound in the ambient medium estimated in Section 14.21 
The relatively small £ value corresponds to a distance from 
the SNR center to the pulsar of s ~ 4 x 10 19 di.5 cue This 
requires the SNR must be elongated along the line of sight 
if J0617 is still within the SNR. For instance, if the SNR is 
a prolate spheroid whose polar axis is parallel to the line of 
sight, then this axis must be a factor of 2 longer than the 
equatorial radius and the length of the chord through the 
SNR along the line of sight to J0617 becomes Zig ~ 7<fi.5 
(cf. Section B~21) . 

In this geometry the north-south jet-like features could 
be true pulsar jets (e . g., sim ilar to the inner jets of the Vela 
pulsar; iPavlov et al.l 20031 1 of ~ 10 17 di. 5 cm (deprojected) 
length. The fact that the southern jet is brighter than the 
northern jet can be explained by the Doppler boosting, 
i.e., the outflow in the southern jet is approaching the ob¬ 
server, which implies that the pulsar is also approaching 
and located in the front part of the SNR. The offset of 
the apparent ring center from the pulsar (see Figure © 
could be explained si milar to that of the Crab inner ring 
dWeisskopf et al.L l2012h ■ e.g., as caused by an intrinsic az¬ 
imuthal asymmetry due to nonuniformity of the ambient 
medium or even a simple geometric displacement caused 
by a non-equatorial ring. 

In this interpretation the observed X-ray PWN beyond 
the ring is comprised of a shocked pulsar wind outflow. 

Its cometary shape is formed by the ram pressure caused 
by the motion of J0617 relative to the ambient medium. 

The misalignment of ~ 50° between the symmetry axis of 
the comet-shaped PWN and the direction toward the SNR 
center could be explained by bulk motion of the SNR mat¬ 
ter in the vicinity of J0617 with a velocity of a few hundred 
km s -1 (e.g., due east or east-southeast). 

We should note that the ellipticity of the J0617 ring 
derived from the image analysis (Section 3.1.2) is subject 
to large systematic uncertainties due to the ring faintness 
and the (apparent) ring-jet intersections. Therefore, we 
cannot exclude the possibility that the true ellipticity is 
much smaller than 30°, i.e., the ring is nearly circular. 

In this case, if we assume that the pulsar’s spin axis is 
aligned with direction of motion, we would have to con¬ 
clude that J0617 could not originate at the center of the 
IC 443 SNR. However, if we allow a strong misalignment 
between the spin axis and t he direction of mo tion (e.g., 
similar to IGR J11014-6103; iPavan et alll2014h . then the 
circular shape of the ring would not impose restrictions on 
the direction of motion and birthplace of the pulsar. Nev¬ 
ertheless, if the jet-like features are indeed the jets, they 
would lie nearly along the line of sight and would be excep¬ 
tionally long. Thus, some ellipticity of the ring is required 
if it is associated with the equatorial pulsar wind. 

The cometary shape of the larger PWN may offer some 
additional guidance in interpreting the observed proper¬ 
ties of J0617. In general, the pulsar wind flows non- 

2(1 Note that the parabola becomes even ‘wider’ if there is a component 


radiatively until it encounters the TS, while the shocked 
pulsar wind, confined between the TS and a contact dis¬ 
continuity (CD) that separates the shocked pulsar wind 
from shock-heated ambient material, emits synchrotron ra¬ 
diation. The shapes of the TS, CD and the outer (forward) 
shock, which separates the shocked ambient gas from un¬ 
shocked ambient gas, depend on the Mach number. 

In the supersonic case, A4 1, the shapes of these 
surfaces ahead of the moving pulsar resemble paraboloids 
whose sky projections look bow-like if the angle between 
the pulsar velocity and the line of sight is large enough. 
The shocked pulsar wind at the bow apex is confined 
in a thin layer, and thus the TS and CD surfaces up¬ 
stream from the central source are close to each other, 
i?Ts(0) ~ Rcd{9) at 9 <C 1, where 9 is the polar angle 
with respect to the symmetry axis of the flow. The TS 
standoff distance can be estimated from the balance be¬ 
tween the bulk wind pressure inside the TS and the ram 
pressure of the ambient matter: 

E/[^cR 2 ts (9 = 0)] ~ pv 2 0 (1) 

(assuming the pulsar wind is isotropic inside the TS). At 
arbitrary 9 the pressure balance defines a surface R ts( 0). 
An analytical solution for Rts(9) was derived by iWilkinl 
(Il996f) in the thin-layer appr oximation, it agree s well with 
numerical simulations (e.g., iBucciantini et all [2005:) for 
small 9 and high M. For 9 <C 1, the Wilkin formula can 
be approximated by the parabola ( y/Ro ) = C(x/Rq) 2 + 1 
with C = 0.3, where Rq = -Rts(O), — y is along the di¬ 
rection of motion, ( x/y ) = tan(9, and the pulsar is at the 
focus x = y = 0. Although the head of the observed 
cometary PWN resembles a parabola, the fit to the X-ray 
contours in Section 13.1.31 gives a substantially larger co¬ 
efficient of the quadratic term, C = 0.54 instead of 0.3. 
The sharper (narrower) head of the observed PWbPl and 
the smooth surface brightness profile ahead of the pulsar 
(i.e., the lack of any structure that could be identified with 
a TS; see Figure 5) provides additional evidence that the 
pulsar motion is not supersonic. 

In the transonic case, A4 ~ 1, the ambient thermal pres¬ 
sure contributes, in addition to the ram pressure, to the 
right hand side of Equation ©• Since c 2 = 7Pthm/p, the 
total pressure can be written as pthm + pv q = pv q(1 + 
7 _1 AI~ 2 ), and the stand-off distance for the CD can be 
estimated as Rcd(O) ~ i?xs(0) ~ E 1 / 2 { / i-Kcpv‘^)~ 1 ^ 2 {\ + 
7 ~ 1 M~ 2 )~ 1 I 2 ~ 10 17 (E/3 x 10 36 ergs _1 ) 1 / 2 plg' /2 (l + 
7 M 2 )- 1 / 2 cm. Inserting plausible parameter values for 
the pulsar and the ambient medium (estimated in Sec¬ 
tions im and 14. 2D . -Rcd(O) ~ 10 17 cm which corresponds 
to 4.5"dj(5 or roughly the observed ring radius. In this 
transonic case, there is no longer a general analytic solu¬ 
tion, but the parabolic approximation for the CD head at 
7AI 2 > 0.7 gives C « 0.3 + 0.2(7.M 2 ) -1 . It becomes close 
to the measured C = 0.54 at yA4 2 « 0.8, or M ~ 0.7 for 
7 = 5/3. 

In the subsonic limit, M < 1, the canonical bow- 
shock model must be abandoned because the thin layer 
assumption, i?xs ~ Rcd, does not apply even in the 
forward direction. Instead, the flow can be treated as 
incompressible so that Laplace’s equation for the veloc- 

of the motion along the line of sight, i < 90°. 
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ity field, with appropriate boundary conditions, correctly 
describes the hydrodynamics. Near the pulsar, the flow 
is again idealized as purely radial, whereas at large dis¬ 
tances the velocity must be uniform and directed oppo¬ 
site the direction of motion of the pulsar. The solution 
is then the sum of a uniform flow and a radial flow from 
the pulsar with speed falling as the square of the radius. 
Since the post-shock radial velocity component must be 
abou t c/3 at the TS outer boundary dKennel fc Coronitl 
11984 but ~ vo <C c/3 at the CD (where the two ve¬ 
locity components cancel), 1?ts(0) -Rcd(O), and the 
CD standoff distance is 7?cd(0) = (c/3vo) 1 ^ 2 Rts(0) — 
18(uo/300kms^ 1 ) _1 / 2 f?Ts(0). We do not clearly see the 
CD in our images, perhaps because it is blurred by insta¬ 
bilities. However, if we assume 1?cd(0) ~ 15"-30"(see Fig¬ 
ures 5 and 6), then the above estimate gives i?xs(0) ~ l"~ 
2"(for vq = 300 km s -1 ), much smaller than the ring ra¬ 
dius, ~ 5". This means the ring may be analogous to the 
Crab’s torus while the TS may be hidden in the wings of 
the pulsar’s PSF (see Figure Note that the relative 
motion cannot be highly subsonic for reasonable estimates 
of vq since c s ~ 400 km/s. 

The overall shape of the PWN also depends on the Mach 
number. In th e supersonic and t r anson ic regimes, numer¬ 
ical solutions (iBucciantini et all l2005h show the Rts(0) 
surface develops a bullet shape with the TS elongated 
in the downstream direction (0 = n) and terminating in 
a Mach disk oriented perpendicular to the direction of 
motion. Simulations show the ratio f?Ts( 7r )/^?Ts(0) ~ 
(1+A4 2 ) 1 / 2 for lower Al, reaching a limit of ~5 for A4 1. 
This amount of elongation is incompatible with observa¬ 
tions of J0617, again suggesting the flow is not supersonic. 
Note that a circular ring is only conceivable if the pul¬ 
sar wind was purely radial producing a limb-brightened 
spherical emission region outside the TS. However, pul¬ 
sar winds are known to be highly aspherical resulting in 
emission morphologies like the Crab torus. As we men¬ 
tioned above, if the ring is a torus, then the spin axis is 
predominantly along the line of sight. 

Finally, in the supersonic and transonic regimes, the ex¬ 
pected shocked pulsar wind sp eeds in the tail ( r oughl y 
along the CD) are 0.3—0.5c (|Bucciantini et all l2005l h 
whereas they are only of order vq ~ 300 km s -1 in 
the subsonic case. Particles in this flow will only ra¬ 
diate at X-ray energies, E ~ 1 keV, for a time r ~ 


200(E/1 keV) _1 / 2 (B/50/rG) _3 / 2 yr due to their syn¬ 
chrotron losse s. Here, B i s the equipartition magnetic 
held defined (iPacholczvkl . 119700 as B = (67rci2(l + 
k)L/(fV )) 2 / 7 where / is the filling factor in a volume V 
producing a synchrotron luminosity L over some spectral 
range, C12 a constant dependent on this spectral range and 
the spectral slope (index), and the ion to electron energy 
ratio k = 0 for a pure pair plasma. From the flux den¬ 
sity in the radio band, ~86 mJy (section [3.2.61) . a volume 
of that region of 10 54 cm 3 assuming a filled prolate el¬ 
lipsoid, and a flat radio spectrum extending to 100 GHz 
(C12 = 7x 10 6 ), the held strength is B ~ 65(sini) 2 / 7 /zG. 
The length of the J0617 nebular tail extends roughly 
1.5' or 2 x 10 18 di.5(sini) -1 cm implying a velocity of 
~500(sin*) -5 / 7 km s _1 for this synchrotron age. This is 
somewhat larger than the subsonic velocity estimate and 
only approaches the transonic velocity estimate as i —> 0°. 

In summary, the shape of the forward region of the J0617 
PWN and the size of its surrounding ring implies a tran¬ 
sonic (or mildly subsonic) net how because of the addi¬ 
tional how collimation that occurs at low A4. The size of 
the circular ring is compatible with the standoff distance 
of a termination shock unless the how is subsonic in which 
case the ring is larger than the TS by a factor of 20 or 
so. In any case, the near circularity of the ring implies 
either emission from a limb-brightened spherical emission 
region or emission from a toroidal region with symmetry 
axis nearly along the line of sight. Neither of these two 
conclusions is satisfactory, however, as the former belies 
the known predominantly azimuthal emission geometry as¬ 
sociated with rotating magnetized neutron stars and the 
latter implies a pulsar orientation unlikely to produce the 
observed shape of the larger cometary PWN. 
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Table 1 

Outlier Regions Spectral Fit Parameters 


Region a 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Area (aremin 2 ) 

4.75 

1.96 

1.53 

3.80 

2.95 

4.23 

5.07 

2.02 

3.03 

Absorbed (Power Law + Variable MeKaL) b 

n H /10 22 (cm 2 ) 

0 49"*" 0,03 
u -^ y -0.05 

n 0^7+0.03 

uo< -0.04 

0 64”*” 0 ' 05 
u -°^-0.12 

0 51+0 0 4 
U.Oi_o 03 

0-48+ 34 

0-48l+ 3 4 

n ctc+0-06 

U-OO-o 06 

n zl &”*” 9,99 
9-48 — 0 05 

0 45~*”°' 03 
u.^o_o 03 

r 

0 ag+°- 33 
o.62 

4 14+1-03 

q 59+O.6O 

o.oz—o 60 

3 91+ 0,51 
-0.30 

o 59+0-08 

O.OZ-i is 

3-50+1° 

9 1 7+0-10 
-0.11 

1 95~*”°' 10 

±.»u_o 09 

2.26+;- 

kT (keV) 

0 69" 1 " 0 ' 02 
u,u -0.02 

0 70+ 0 03 

u - * u _o.oi 

0 79+ 0 05 
u ' ' y -0.16 

0 69"*"°' 02 
u,u -0.02 

0.68+°; 33 

0 66"*" 0 02 
u.uu_o 02 

0.681+1 

O 70+0.06 

U '' 0.06 

0.68++ 

Xpow/10 -4 

1 41+0.99 
0.92 

0-04+;“ 

1 30+ 0 - 73 
±.ou_o 96 

2 fifi+ 1-51 
z -OO_0.96 

0-85++ 

1 71 +0-63 

L -' 1 -0.95 

3.341+1 

2-061+1 

2 45 +0 - 42 

z, ^°—0.51 

Xthm/ 10 -4 

19 If)”*” 5 ' 11 
1Z. 1U_ 2 48 

8-20++ 

q o£!+3.52 

o.zu_ 2 28 

8 34+ 3 - 38 

0.o^_3 71 

3 01"*” 1 ' 64 
°' u± -1.16 

Q C7+3.01 

0m01 -0.23 

9 qq+ 1, 09 
^•oo_0.35 

0-87++ 

4 i0+ 803 
^•J-O-o.51 

/pow/ 10- 13 ( c ) 

1.60 

0.30 

0.88 

1.75 

0.70 

1.42 

7.58 

6.23 

5.43 

/thm /10 13 

8.59 

6.74 

1.57 

6.43 

2.68 

4.31 

3.18 

1.04 

5.32 

X 2 /dof 

356.6/238 

251.0/180 

176.6/161 

346.4/222 

314.2/188 

359.7/220 

507.0/278 

307.9/259 

381.9/262 

Absorbed (Variable MeKaL + Variable MeKaL) d 

n H /l0 22 (cm 2 ) 

q 47+0.04 
* —0.05 

n qq+0.03 
u.oo—o 03 

0 53+0-15 
u ' oo -0.10 

q 47+0.07 
U '^'-0.05 

0 45+ 0 - 08 
u>q:O -0.09 

0-48+06 

0 50+ 005 
U-OD_o 06 

0 49+O O 9 
u ‘^ y —0.08 

0 46"*”°' 04 
u.^O-o 04 

kT i (keV) 

0 6O"*" 0 04 

U.UU-O 04 

0 65“*" 9 - 93 
u -°°_0.03 

0 5Q+ 0 - 07 
U.OU-o 07 

0-60+° 4 

0 5fi+ 0 - 08 

U.OD-o 08 

0-551+1 

0 64"*" 0 02 
u ' u -0.02 

0 64"*"°' 05 
u,u -0.05 

0 07+0-02 
u,u ' -0.02 

fc+2 (keV) 

n 70+0.02 

U ' ‘°-0.02 

0 74+°' 15 
u ' ‘ -0.02 

1 02”*”°' 26 
±.uz_o 09 

0-80i 

0 77+0.11 

U '' '-0.05 

0.821+1 

q cc+0.36 
O.OO-o 30 

4 72+ 0 - 78 

** '-0.50 

q 75+1-00 
°-'°-0.51 

Ab/io- 4 

3 19”*” 3 ' 18 
0 '- Lt7 —0.29 

9 £10+0-32 
^•O»_ 0 .57 

4-27+;“ 

7 44+ 6 - 40 
* '”—4.07 

1 79+1.57 

A - ' z -1.08 

3-601+1 

6.81I+® 

2-98+°7 

7.39+ 43 5 

X2/IO- 4 

16.86+;“ 

9 13+ 2,60 

».xo_ 7 19 


12 19”*” 4 ' 94 

iz.iy_7 16 

6-!8 1+7 

6-601+ 4 

7-101+1 

5.68++ 

4-01++ 

/i/10~ 13 

2.52 

2.39 

0.70 

3.37 

0.88 

2.32 

3.43 

1.22 

5.97 

/2/IO- 13 

7.41 

4.37 

1.70 

4.73 

2.43 

3.28 

7.09 

5.87 

4.79 

\ 2 /dof 

323.1/238 

239.7/180 

150.1/161 

327.3/222 

305.3/188 

317.5/220 

414.5/278 

280.1/259 

374.4/262 


a Numbers correspond to the numbered regions shown in Figured] 

b Abundances tied to solar ratios for the groups CNO, Ne and Mg, Si and S, and Fe and Ni. 
c All fluxes in units of erg cm -2 s -1 . 

d Same as (b) and abundances of second MeKaL model tied to those of first model. 


Table 2 

Spectral Fit Parameters & Derived Quantities for CXOU J061705.3+222127 


Model 

7lH 

10 22 cm -2 

rj~< OO 

eV 

R°° 

km 

T OO 

-^bol 

10 32 erg s 1 

r 

Lt 

10 31 erg s 1 

X 2 /dof 

NSA 

BB 

0 63 +u,U2 

U.OO-0.03 

0 61 + 013 
u - ol -0.09 

58.41+ 

1.32.01+ 

13.05 

1 co+0.09 
1 - oo -0.07 

2.6±0.1 

1 o +0 - 2 

1 Q4+U-31 

r.OT_ 0 .35 

2 00+°- 30 
Z.UU_o.20 

1 o +i - 4 

4 - u -0.5 

1 4 + 0.6 
4 -^-0.5 

42.2/42 

41.7/41 
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Table 3 

Upstream Spectral Fit Parameters 


Region 0 1 2 3 4 5 6 


Absorbed (Power Law + Variable MeKaL) b 


nn/10 22 (cm 2 ) 

o-6il8:8l 

0 . 40+04 

n 4Q+ 0 - 06 
u -^ u - 0.05 

n 4o+ 0 05 
u -^ u - 0.05 

0 40 + 0.08 
u -^ u - 0.06 

fl 97 + 0.06 
u -°' - 0.05 

r 

2-08++ 

9 1 Q + 0.23 
z - iO - 0.24 

2 1Q+ 0 - 27 
z - iy - 0.28 

2 5Q+ 0 - 39 
Z . O^_038 

2 41 +°- 44 
Z '^ 1 - 0.47 

2 4Q+ 0 - 41 
Z '^ c > - 0.42 

kT (keV) 

0 9Q+ 0 - 10 
u . oy_ 0> 05 

0-70+88 

0 67 +0 - 07 
u -°'- 0.10 

0 72+ 0 - 12 
u -' z - 0.04 

0 68 +0,07 
u . oo _ o.08 

0 71+ 0 - 10 
u - 1 i - 0.09 

-Kpow/10 -5 

6 50+ 115 
O . OU - 0.98 

3.oi±8:®T 

2 4^+ 0 - 82 
Z '^°- 0.66 

2 9Q+ 1 - 01 
0.78 

1 crq + 0.89 

1 62 +0,79 
i * oz - 0.61 

/' i hm/IO " 5 

7 i o + 3.76 

1 - lz - 2.53 

9 yq + 0.93 

1 ^- 0.79 

2 Q1 +0, I 9 

z - yi - 0.70 

9 7 c ;+ 0-71 
z - ‘ °- 0.64 

o 09 + 0.81 
°- oz — 0.71 

Z - UO - 0.59 

/pow/ 10 -13 (erg cm” 2 s” 1 ) 

1.60 

0.78 

0.60 

0.41 

0.32 

0.32 

/thm/ 10” 13 (erg cm” 2 s” 1 ) 

0.26 

0.31 

0.23 

0.22 

0.26 

0.23 

X 2 /dof 

135.5/146 

120.1/116 

111.7/ 99 

97.8/ 90 

58.6/ 84 

66.5/ 80 


Absorbed (Variable MeKaL + Variable MeKaL) c 


nn/10 22 (cm 2 ) 

0.63t8+ 9 

n 71+0.10 
u - 1 -*--0.13 

n 7Q+ 0 - 14 
u -' y -0.11 

n 07+O.O6 
u -°‘ -0.04 

0.40+og 

0 9f)+°- 07 
u - oo -0.05 

fcTi (keV) 

0.38+05 

0 42+ 0 - 07 

U-OZ —0.03 

0 29 +0,04 

u - zy -0.05 

q 71+0.05 
u *' 4 —0.03 

o.67i8;8? 

n 7Q+0.07 

u -' u -0.09 

kT 2 (keV) 

4-i6+:?8 

9 c;c;+0.51 
Z ' OO -0.39 

2 40+°- 56 
z -^ u -0.38 

q q7+1.43 
-1.09 

4.37+8? 

9 09+ 3 - 23 
O.UO-i 23 

Ri/10” 4 

0.90+j 4 

1 91 +0-81 
1 -° ± —0.59 

1 6^+ 2 - 24 
i.OO—0.71 

f) qq+0.07 
O.OO-0.04 

o.39i8;8? 

0 90 + °- 07 
U.OU_ 0 .o7 

K2/IO- 4 

1 86 +011 
-L-OO-o.io 

n 07+0.08 
u -°‘-o.io 

0 69 +0,08 
u.oy_ 0 .o7 

n qi + 0.O8 

u - Oi -0.05 

0 23 +0 - 07 
u - zo -0.06 

n 24+ 0 - 07 
u - z ^-0.06 

/1/IO” 13 (erg cm” 2 s” 1 ) 

0.31 

0.30 

0.24 

0.28 

0.31 

0.26 

/2/IO” 13 (erg cm” 2 s” 1 ) 

1.54 

0.74 

0.55 

0.36 

0.29 

0.29 

X 2 /dof 

144.1/146 

134.0/116 

115.2/ 99 

107.9/ 90 

65.5/ 84 

72.3/ 80 


a Regions upstream (southwest) of J0617 (see Figure [8]), numbered by increasing distance from J0617. 
b Abundances fixed to those of Outlier region 6 (see Table [Q 

c Same as (b) and abundances of second MeKaL model tied to those of first model. 
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Table 4 

Pulsar Wind Nebula Spectral Fit Parameters 


Region 3, 

1 

2 

3 

4 

5 

6 

7 

Area (arcmin 2 ) 

0.048 

0.091 

0.157 

0.266 

0.365 

0.531 

0.722 

Absorbed Power Law 

uh/IO 22 (cm 2 ) 

0 %+°- 04 
U -°°-0.04 

f) crq+0.03 
U ’°°-0.03 

0 ^fi+O- 04 
U -°°-0.03 

0 ^9+ 0 - 03 
u ‘ oz -0.03 

0 +0-03 

u - ol -0.03 

n 47+0.03 
-0.03 

q 41+0.02 
u '^ ± -0.02 

r 

1 ^n + 006 

10U -0.03 

1 54+ 006 
r- J ^-0.05 

1 69 +006 

1 0.03 

i-raigS 

1 Q1+0-06 

1 -^ 1 -0.06 

2 06 +0,06 
z - uu -0.06 

9 9Q+0.07 
z ' zo -0.07 

ATpow/ 10 -4 

1 qr+0.10 

i.oo_oo9 

1 ^ 4 + 0.11 
i -°^-o.io 

1 7Q+ 013 
* y -0.11 

1-8418:18 

1 Q«+°- 14 
i-yo_ 0 .i3 

2 02 +0 ' 14 
z - uz -0.13 

2.0818:1 

/pow/ 10 -13 (erg cm " 2 s -1 ) 

6.84 

7.56 

7.20 

6.72 

6.20 

5.42 

4.61 

/pow/iO " 13 (erg cm " 2 s -1 ) 

9.14 

10.11 

10.21 

9.71 

9.40 

8.61 

7.78 

X 2 /dof 

289.2/293 

336.7/302 

278.9/295 

326.5/285 

294.4/282 

353.2/274 

375.6/254 

Absorbed (Power Law + Variable MeKaL) b 

n h/ 10 22 (cm 2 ) 

0 f)4 +0,35 

U '°^-0.07 

0 67+ 0 - 25 
u - 0 ' -0.16 

0.75/°+ 

n 77+0.18 
u -' * —0.15 

0 07+ 0 - 18 
u - u ' -o.io 

0 60 +0 ' 15 
w.UU-o.09 

0 47+0- 09 
-0.07 

r 

1 cc+0.09 
-L-OO-o.04 

1 crq+0.03 

±.oo_o.o3 

1 77+O.I4 
1 *‘ ' —0.13 

1 «Q +°- 05 

i-oy_ 0 .il 

1 qc;+ 0- 12 

± ' yo -0.09 

2 03+ 012 
Z.UO-0.09 

1 07+O.IO 

1 -0.09 

kT (keV) 

0 1 1 + 015 
u<li -0.03 

U.OU-o.28 

0 %+ 5 - 12 
u -°°-o.io 

0 qo+0-09 
U.OZ-o 04 

0 OQ+ 0 - 04 
u - zy -0.05 

n 97+0.03 
u< °‘ -0.06 

0 67 +0 ' n 
U.U/-006 

A"pow/10 ~ 4 

1 45+0.31 
l-W-0.13 

1 4Q+°-47 
i -^ y -0.47 

2 03+°- 45 

Z.UO-037 

9 on+0.36 
z.zU-0.33 

2.14l° : 3 ° 

i-98l8:l 

1 41+0.47 

4 •^ 4 — 0.16 

ATthm/10 -4 

6Q q4+ 3837 - 17 

1 02 +6 ’ 70 
±.UZ _1 02 

0 ^ c i + 4 - 52 

U.00-0.27 

0 . 79 ++ 

1.891+7 

0-6918:1 

0 61 +1 ' 14 
U - U1 -0.12 

/pow/ 10 -13 (erg cm -2 s -1 ) 

6.81 

7.56 

6.99 

6.44 

5.98 

5.20 

4.10 

/thm/10 -13 (erg cm -2 s -1 ) 

0.00 

0.00 

0.16 

0.24 

0.21 

0.28 

0.70 

/pow/ 10 -13 (erg cm -2 s -1 ) 

9.40 

10.10 

10.85 

10.60 

9.81 

8.62 

6.39 

/thm/10 -13 (erg cm -2 s -1 ) 

0.80 

0.08 

2.03 

3.57 

2.46 

1.62 

1.98 

X 2 /dof 

284.8/287 

336.7/301 c 

267.2/289 

298.0/279 

258.8/276 

276.3/268 

215.5/248 

Absorbed (Variable MeKaL + Variable MeKaL) d 

n h/ 10 22 (cm 2 ) 

0-54l°;°o5 

f) cq+0.07 

U.OO_ 0< 06 

0 - 65 /°+ 

0-5118+5 

0-6918+ 

n 44+0- 10 
0.04 

n 41 + 0.07 
U '^ ± -0.06 

fcTi (keV) 

0.08™ 2 

1 f ) S +78 ‘ 87 

l.UO-o.79 

0 . 63 /°;- 

o.46l8;18 

0 . 7118:85 

o-4il8:8i 

0 (l c i+ 0 - 06 
u - oo -0.06 

kT 2 (keV) 

14-70++ 

23 60+ 56 - 30 
zo.ou_ 19 58 

9-n±?;iS 

7.86/1:88 

7-981+8 

4.981+11 

r q'7+1-23 
001 -0.81 

AV/10 -4 

one 70+342.30 
OUO. ( O _ 305 73 

1 6 7+ 2 - 49 

1-59 ±lf 7 

1 oq+5.53 
4 '°^-1.32 

O qr + 1.72 
0.09-1 66 

9 00+3.76 
Z.ZZ-1 16 

2-691+? 

k 2 / 10 -4 

c 70+O.4O 
* °-0.95 

cj Q^+0-34 
° >JO -5.95 

5.67++ 

3 89 +0,44 
o.oy_ 0 .36 

4 44+ 0 9 2 
^•^-0.87 

e 9C+O.37 
°- zo -0.52 

q 01 +0.35 
°’ O1 -0.40 

/i/10 -13 (erg cm -2 s -1 ) 

0.38 

0.91 

1.70 

0.75 

4.25 

1.02 

2.18 

/ 2 /IO - 13 (erg cm -2 s -1 ) 

6.35 

6.62 

5.41 

5.84 

1.97 

4.30 

2.54 

/COT/iQ- 13 ( erg cm -2 s -i) 

0.38 

0.91 

1.70 

0.75 

4.25 

1.02 

2.18 

/™ r /10 -13 (erg cm -2 s -1 ) 

8.74 

9.11 

9.44 

8.63 

7.65 

6.96 

5.32 

X' 2 / dof 

280.0/287 

337.3/296 

266.0/289 

301.4/279 

254.2/276 

277.5/268 

212.4/248 


a Regions ordered outward from J0617 numbering 1 through 7, respectively. 
b Abundances tied to solar ratios for the groups CNO, Ne and Mg, Si and S, and Fe and Ni. 
c Abundances frozen to best-fit values before errors on free parameters estimated. 
d Same as (b) and abundances of second MeKaL model tied to those of first model. 
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Table 5 

Pulsar Wind Nebula X-Ray & Radio Spectral Fit Parameters 


Region 0 

1 

2 

3 

4 

5 

Area (arcmin 2 ) 

0.05 

0.10 

0.14 

0.22 

0.34 

7th/10 22 (cm 2 ) 

0 4f)+ u - u3 
u -^ u -0.03 

0 ^+ U - U4 
U.OO-0.03 

0 ^7+ u - U3 
u -°‘ -0.03 

0 52 +u ' u:;} 

u - oz -0.03 

f) +U.08 

u - J1 -0.03 

r 

1 47+0.05 
-0.05 

1 61+ 0 - 06 
± - ul —0.06 

-1 70+0.06 
- 1 -* ‘ °-0.06 

1 01 +0.06 
1 -° i -0.06 

1 QQ+0.06 

±.^o_o 06 

A'pow/lO -4 

1 OO+0.08 

1 0.08 

1 5«+ 0 - 12 
i - oo -0.11 

1 7 , '7+0.13 
' —0.12 

1 54+ 012 

-L-O^-o.n 

1 QQ+O-i-3 
1.00-012 

/pow/10” 13 (erg cm” 2 s” 1 ) 

7.67 

7.02 

6.74 

5.76 

5.73 

/pow/10' 13 (erg cm” 2 s” 1 ) 

9.65 

9.61 

9.73 

8.42 

8.76 

x 2 /dof 

299.6/305 

271.8/289 

274.1/285 

285.5/269 

370.3/268 

/330MHz ( m Jy) 

5.0±0.9 

11.0T3.1 

11.5±5.6 

10.5T8.4 

18.0T13.0 

v (GHz) 

0 02 +0 16 
u.uz-o 02 

l.nttl 

cyj -\-66 
z ‘ -20 

128/gg 8 

826/479 


Elliptical annuli numbered from inner to outer 1 to 5 







